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Abstract. Using the modified factorization method employed by Mielnik for the 
harmonic oscillator, we show that isospectral structures associated with a second 
order operator H, can always be constructed whenever H could be factored, or exist 
ladder operators for its eigenfunctions. Three examples are shown, and properties like 
completeness and integrability are discused for the general case. 



1. Introduction. 

In the context of one or tree-dimensional Quantum Mechanics, the principal problem 
is to solve the eigenvalue problems associated with the Schrodinger equation when the 
particle is under the influence of a general potential. The two most commonly used exact 
methods to determinate the eigenfunctions and the eigenvalues (once we have separated 
the differential equation) are the method of orthogonal functions and the factorization 
method. Nevertheless, the number of potentials which have exact solution by these ways 
is quite small. 

In recent years, Abraham and Moses [Q] showed one algorithm that allows us to 
generate exactly solvable potentials from another potential which has exact solutions, 
both with the same eigenvalues, but adding or subtracting a finite number of them 
through the systematic use of the Gelfand-Levitan equation [§]. 

On the other hand, Mielnik || observed that for the harmonic oscillator the same 
isospectral class of potentials that Abraham and Moses obtained, may be easily built 
using the ladder operators. On the same way, Nieto |4j showed that into supersymmetric 
Quantum Mechanics, the natural factorization that of the bosonic or fermionic one- 
dimensional Hamiltonian could produce isospectral Hamiltonians too. Nevertheless, the 
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use of the factorization scheme to produce new classes of solvable potentials from other 
one is not restricted to the two kinds before, as we will show. 

The present work has been divided into the following sections, in section 2, we 
resume the procedure used by Mielnik in order to find an isospectral family to the 
harmonic oscillator. Along section 3, we show how this property is not exclusive of that 
case, but all the second order differential operators which may be factored by ladder 
operators have it. Moreover, this property exists always that a second order operator 
can be factored in a more general way, not necessarily by ladder operators or in a kind of 
factorization like the supersymmetry one. Section 4 contains three different examples of 
this techniques, the free particle in one and three dimensions and the isotropic oscillator. 
Finally, our conclusions are summarized in section 5. 



2. The oscillator's isospectral family 

The first family of isospectral operators built in a different way from that of Abraham 
and Moses [U , was obtained for the one-dimensional harmonic oscillator, using a variant 
of the factorization method proposed by Mielnik [[J . In this case, the Hamiltonian 

Id 2 1 2 , . 

H = ^ + -x 2 2.1 

2 dx 2 2 V ; 

whose eigenvalues are E n = n + ^, is factored by the creation and annihilation operators 



in such a way that H = a^a + \ = aa) — |. The method considers the more general 
operators 

b =T2^ + m )' ^7fK +/J( 4 (2 ' 3) 

which must be such that H = bb* — |. 

Moreover j3{x) satisfy the Ricatti equation j3'+j3 2 = x 2 + l, which have the particular 
solution p = x. Then the general solution can be obtained putting (3 = x + <f>(x). This 
yields 

<P(x) = "r — | — rr- (2.4) 
where A is a parameter. Now, we write down the inverse product tfb 

H x = b^b+ l - = - l -+ X -x 2 - <j/(x) = H — <j)'. (2.5) 
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Hence, H\ is a certain new family of Hamiltonians. Moreover, if we define ip n = b'(p n -i, 
where (p n are the eigenfunctions of H, thus, ip n are the eigenfunctions of H\, with 
eigenvalues E n = n + ~ 

ff A V„ = 6 t ( J ff+l)^_i = (n+i)Vn. (2.6) 

Therefore, if and _£/a are isospectral, except for one missing element, the ground state 
of H\, which should be orthogonal to all of ip n 

WoM = 0^6^0 = 0. (2.7) 

We will call to {H, <p n } and {H x ,ip n } isospectral structures (IS). Some other usual 
properties of the harmonic oscillator are inherited for its IS, like the coherent states, as 
has been shown before ||. As it can be seen immediately, this is not the only case for 
which there must exist IS |J. The existence of IS for a second order differential operator 
H is related with the existence of a factorization scheme, as we shall show below. 



3. IS by the Factorization Scheme. 

We begin by establishing the immediate generalization of the Mielnik's results. Consider 
the second order differential operator 

H = P(x)^ + Q{x)^- + R(x) (3.1) 
dar ax 

and its eigenvalue equation 

Hip n = \ n Lp n . (3.2) 

By simplicity, we will supose that n is an discrete index, and that A n has a minimum 
value. We assume that there exist the first order ladder operators, then we may write 
them as 

where A^ip n = c£(p n ±i, with constants, and they are such a way that factored to H 
into 

H = A + n _,A- + k\ (3.4) 

= A n+1 A^ + k 2 n , (3.5) 

the constants k\ and k\ are related with A n by AA n = A n+ i — X n = — k\. Notice 
that of equations ( |3.3|) , (3.4) and (3.5) a n , (3 n , 7 n and 5 n satisfy 

— Pn-lOn — Pi 

Pn-lln - Pn-l&n ~ a n-lK = Q, (3-6) 
a n -lln + Pn-ll' n + K = R, 

a n _i7„ - b~ n ot n _ x + k 2 n = R. 



Next, consider the more general operators 



a { n ] = &n{x) + Pn-j-, «„ 5 = 7„(x) - S n -^-. (3.7) 
dx dx 

obtained by the changes a n — > a n = a n + T] n (x) and 7„ — > % = 7„ + u n (x) in Q3.3|), with 
7? n and z/ n functions. Then, we will have the two following cases 
Case I. Assuming that H is factored by into the form 

H = 0S0H + kl (3.8) 

and using the equations ( |3.7D and (3.6) we obtain for rj n and v n 

Ai-i 



Vn-l 



and the Ricatti equation for v n 



'A 



ra-1 



'A 



'n-l 



whose general solution is 



A+ / dx 



X„ = dx 



®-n— 1 _|_ Tra 
An-l n 



(3.9) 



(3.10) 



(3.11) 



Hence, the factorization of H into the form ( |3.8| ) is not unique. Moreover, building the 
inverted product an+\ a n^> we may find the new one-parametric family of operators 



(3.12) 



which as we can immediately see, has identical spectra to that of H, even though their 
eigenvectors are different. If we define 



then 



H*ip n = a n+ \ (H - AA n ) (p n+ i = \ n ifj n . 



(3.13) 



Therefore, ip n are the eigenvectors of and one state could be added to the spectra 
A n - 

Note that in the more general case defines a different operator for different n. 
Hence, to be precise, H and share only one element of their spectra: A n . In fact, 
the functions {ip n } are not orthogonal in general, because equation ( |3.13j ) gives 

(A n - A m ) (ip m , ijj n ) = (ip m , A mn ip n ) , 

with A mn = — H^, a non vanishing scalar function. However, when the index n is 
different to that defined by equation ( |3.2| ) , like the index / for the angular momentum into 



the Schrodinger radial equation with central potentials, the dependence of in n may 
be seen as a result of the action of one potential, which is coupled in a different way for 
each state, like the effective central potential, with the orthogonality and completeness 
established by the Sturm-Liouville theory M. 

On the other hand, when a!~^ = a^_ 1; then will have the integrability 

properties of {f n }, except for a change on the normalization 

O0n, Ipn) = (An ~ & n ) (<fn, <Pn) ■ 

It is clear that if the ladder operators do not depend on the index n, all the properties 
of {(fin} must be inherited by {ipn}, and probably there is another eigenstate of H x 
obtained as the solution to 

{xM = Q -> a (+) X = 0. 

Case II. The above results may be obtained again, when factored H into 

H = a^ ia W+k 2 n . (3.14) 

In such a case, we find v n+ i(3 n = S n+ ir] n and the Ricatti equation for r) n 



( <Wl \ 

7 " +I + AT 



. / n+ i \ 2 / n 



whose general solution is non trivial, as in the case I. Hence, we may define the new 
family of operators, like ( |3.12| ) 

= a^XaP + kl = H + {f3 n ^v' n + 5 n ^_i} , (3.15) 

whose eigenfunctions 

$ n = aZli^Pn-i (3.16) 

have the same eigenvalues that of ip n . But now, the spectra has not the ground state. 

Note that the comments for the case I still hold in the present one. 

Once we know the existence of IS for H, we must emphasize that the constriction on 
to be ladder operators and factored to H at the same time, is not entirely necessary, 

IS can be built if just one of them holds as we will show below. 

Let be ladder operators, in the more general case they do not factorize to H ||, 

instead, 

D n = A~ +1 A+ and E n = At_ x A- n (3.17) 

define two different operators, which obviously share eigenfunctions with H, but not 
eigenvalues. Even now, the algorithm discussed above can be reproduced to build IS 
associated with {D n , ip n } and {E n , <p n }, but it is important to note that we do not have 
IS for if in a natural way. Nevertheless, when A^ and H are related by 

(D,: ^ f(x)(H-X n )+( dn ), (3.18) 



E n / \ e 



n 
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d n , e n being constants, it is an easy task to find IS for H, as we shall note in the examples 
below. 

Now, if we assume that H is factorizable, not necessarily by ladder operators, then 
there are two first order operators B n and C n which satisfy 

H = B n C n + k n (3.19) 

where k n is a constant. In this case, it is always possible to follow the usual procedure in 
order to find IS for H. Moreover, as F(„) = C n B n is in general a different operator to H, 
IS may be built from it in a direct way, whenever F^C n = C n (B n C n ) = C n H^ n ). Hence, 
£ n = C n ip n are eigenfunctions of F n with eigenvalue A n . Thus, we have established that 
is sufficient with factored one second order operator H , to assure that one IS {H x ,£ n } 
associated with {H, ip n } must exist, even when the spectra is continuos. 

In the case when H does not depend in n, neither B n , C n and k n do, and then the 
family H x will be defined in a proper way. Moreover, in such a case, the completeness 
and orthogonality of {£„} must be insured by the properties of the eigenfunctions of H, 
or by the Sturm-Liouville theory. 

Finally, the present approximation is not the more general. It is possible to 
generalize the four functions a n , f3 n , 7„ and 5 n at the same time, but this scheme 
modifies H in something more than one scalar function, in such a way that a new 
family of isospectral operators to H can be found, but the new operators might have 
not a direct physical interpretation, and just have mathematical sense. These kind of 
IS will not be considered here. 



4. Examples. 

We shall illustrate the above techniques for several examples: the free particle in one 
and three dimensions and the isotropic oscillator. 

4-1. The Free Particle in one dimension. 

We first consider the most simple example: a free particle, that is, one for which 
V(x) =0. The Schrodinger equation is 

d 2 

H(p k = --^^k = k 2 (p k (4.1) 

where k 2 = 2mE/h 2 . Obiously, Lp k are sin(A;a;) or cos (for). H may be factored into 
H = AA* = A^A, with A=£. 

Next, consider the generalized operator 

a = — + u(x). (4.2) 



Thus, if H = aa^ this leads to 

*«) = (4.3) 
Hence, we may propose the new operator 

H* = a^a = -^- 2+ V x (x), (4.4) 

where 

v*(«) s a/<*) = (4.5) 

The above potential has one singularity in xq = —A and vanishes for x — > ±oo. The 
eigenfunctions of if A are 

Xk = aV* ^ A Xfc = ^Xfc- (4.6) 
These eigenfunctions are orthogonal 

(Xfc, Xfe') = (oaVfc, = & 2 (^fc, ( 4 - 7 ) 

Moreover, as ax = has the only solution x = z/(x) and Xo( x ) = v { x )-> then, {xk} is 
a complete set, and like {(fk}, they are an external basis for the Hilbert space. Other 
IS can be found if we assume now that H = a^a, but it is essentially of the same kind, 
because in this case we will obtain v — (A — x)~ l . 

Unlike the algorithm of Abraham and Moses [0 , when we consider the particle in a 
box, the boundary conditions which satisfy the wave functions <p, are not inherited by 
the eigenvectors of the new Hamiltonians. Therefore, the spectra of one particle in a 
box with potential V\(x), must be found by making that the general solution satisfies 
the boundary conditions of the box, which give us different quantization conditions on 
k to that of the box without potential, and then the new spectra will not be equal to 
the particle in a box. 

The present example may be extended straightforward to the three dimensional 
problem. However it is more illustrative the use of the ladder operators along one of the 
indices which does not denote the energy level. 

4-2. The Free Particle in three dimensions. 

The Schrodinger equation for the free particle in three dimensions takes the form 

h 2 

= V 2 V> = Eip. (4.8) 

2m 



Equation ( |4.8|) can be separated in spherical coordinates, introducing ip(r,9,4>) = 
Ri(kr)Yi m (9,<p). The radial equation becomes the spherical Bessel equation 

I A (f*m) + U - <«±V) R M =0, (4.9) 



p 2 dp \ dp I \ p 



where k 2 = 2mE/h 2 , and p — kr, it has as general solution one linear combination from 
the spherical Bessel functions ji(p) and ni(p) 0. Although the radial equation defines 
the energy, it may be written down as an eigenvalue equation for I, hence, the procedure 
showed in the section |3| can be applied here. Consider the recurrence relations [7] 

21 + 1 



fi-i + fi+i 



x 



-ft and lfi- 1 -(l + l)fi +1 = (2l + l)fl 



(4.10) 



satisfied by ji and ni. Using the last equation we may build the ladder operators 
A] 



— and A, 

P dp 



l±l + iL SU ch that Affi = f l±1 (4.11) 
p dp 



which factored Hi into Hi = Af_ 1 A l = A l+1 Af. Hence, the IS associated with Hi will 
be 

I.- The new family of hamiltonians H Xl = Hi + V Xl (p), with 



V l x ^(p) = 2^ 1 (p) and ^(p)=p 21 



dp 



Ai- 



P 



21+1 



(21 + 1) 



whose eigenfunctions are 
*i(p) = 



d l + l 
dp p 



Ri{p)- 



(4.12) 



(4.13) 



As can be seen, if Ai < 0, then neither u Xl and V Xl have singularities, and therefore 
neither ip Xl has. Moreover, as v Xl — > when r — > oo, it assures the integrability of the 
eigenfunctions obtained from ji. 

II.- The other one parametric family is H X2 = Hi — V X2 (p), where 

n -l 



V l x *{p)=2-vf{p) and v\\p)=p- 21 
dp 



A 2 + 



P 



-2Z+1 



-2Z + 1 



and in this case the eigenfunctions are 



i 



dp p 



Ri(p)- 



(4.14) 



(4.15) 



Now, as Vq 2 = (p + A2) 1 , it has a singularity in p = — A2 while i^ A2 has one in p = if 
A2 < 0. For p << 1 the limiting values for vf 2 have the form 

21 - 1 



and therefore the respective limiting form of ip X +i is 



^1+1 



d I 21 + 1 
dp p p 



Ri 



-A l Ri = —R1-1, 



(4.16) 



(4.17) 
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moreover, considering the limiting values of ji |7j 

*W ~ prri^ «•»> 

we note that has not singularities for Z 7^ and A2 < 0, while i/j^ 2 has one in p = 
like 1/p, but it do not appear if we use R = n instead j . On the other hand, as 
vanishes in the limit r — > 00, then the eigenvectors always go to zero at infinity. 



4-3. The Isotropic Oscillator. 



The isotropic three-dimensional harmonic oscillator has the characteristic that when the 
Schodinger equation is separated in spherical coordinates, the eigenvalues of the energy 
are defined in function of two indices, the angular momentum I and that of the radial 
equation, in such a way that IS may be build by ladder operators associated to each 
index. The radial equation for this problem takes the form || 

d 2 



HiUni 



2 Kl + 1) 
+ r 2 + ^ 



U n l(r) = SnlUnl 



with u n i(r) = rR n i(r), e n \ = {An + 21 + 3). The solution to ( |4.19| ) is 

u nl (p)=p {l+1)/2 e-^L l n + Hp) 

where L^ +1//2 (p) are the associated Laguerre polynomials, and p = r 2 
recurrence relations 0, [l(J 

(n + l)L k n+1 = (2n + k + l — p)L k n - (n + k)L k n _ x 

= nL k n - (n + k)L\ 



n d T k 

rfc-l 



7 rt-l 



j k j k 
— L n-1 

■k 



pL™ =(n + k)L k n _, - (n - p)L^ 

we may find the ladder operators over n 



and over k 



A: 



Bl 



n — p 



k + p- 



d 



(n + fc + 1 - p) + p 



dp 



5+ = fi^ 



1 



P- 



(4.19) 

(4.20) 
From the 



(4.21) 



(4.22) 



dp K r dp 

for the L k n functions, which were obtained by J. Morales et al. [11] in a different way. 
Using the last operators is an easy task to build the ladder operators over n or I for u n [. 
Note that we can write 



Uni(p) = fi(p)L l n 2 ; fi 



fl+l)/2 e -p/2_ 



(4.23) 
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Thus, the ladder operators on Z become af = fi±\ B^ +l j 2 f l 1 , which in function of 



may be written down as 



+ r + — 
r dr 



Z + l d 

h r — 

r dr 



(4.24) 



Hence, [a^Y = a^_ 1 , and then the integrability of the new eigenvectors is assured 
from the integrability of {u n i}. Moreover, Hi is factored into H = ajjtiaf — (21 — 1) = 
a T+i a t ~ + 3). As above, there are two class of IS associated with Hi 

I.- The family of Hamiltonians H* 1 = Hi + VJ Al (r), where 



and the eigenfunctions of H l 1 

d 



, +r + - + z/ z Al 
dr r 



[Ai-^xVdx] 



(4.25) 



(4.26) 



with eigenvalues given by e n i. It is important to note that if Ai < then 1 — > when 
r — > 0, oo; and it has not singularities. After some algebra, it can be show that the new 
eigenfunctions have not singularities, and they go to zero when r — > oo. 



II.- The other isospectral class is formed by the Hamiltonians H^ 2 
and its eigenfunctions 



Hi + V, 



*2, 



d Z + l A2 

dr r 



(4.27) 



where now, 



V?* 



-2- 



dr 



,A 2 . 



A 2 



[A2 — J r °° £ 2 'e x dx] 



In this case, if A2 < and Z 7^ 0, then z/ A2 and V* 2 — > when r 



rA2 



singularity in r = like 1/r while 2 diverges like 1/r 2 . Again, when we analize v{ 



(4.28) 

00, but v^ 2 has a 

A 2 



around r = we find 



and then 



"i 2 ~ "2 



r + 



, U n l 
Unl-1 H 

r 



(4.29) 



(4.30) 



but because ( |4.20| ) the eigenfunctions have not divergences. For Z = one can easily 
check that if A2 < or A2 > Uq 2 has not singularities and so neither Vq 2 and ipno 

have. 
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It is important to note, that both cases above are very different to that reported by 



Fernandez [ 12] and Dongpei [13 



An analogous procedure may be followed to build the ladder operators on n for u n i, 
though we should note that when n — > n + 1, the energy level is moved two levels, 
therefore, we must consider two kinds of IS, depending on what part of the spectra we 
are moving. In order to simplify, we shall consider only the case / = 0. So, the ladder 
operators on n may be written as 



r 



1 

n H 

2 2 



r d 
2dr 



r 2 r d 

71+1 1 

2 2dr 



(4.31) 



then (a fl )t = a^-ij ^ assures the integrability of the eigenvectors associated with the 
new structures. Note that in this case, H is not factored by a J, instead we have 

~ 2 - 3^ 



a n+l a n 



H - e n0 \ +(n + l)(n + -„ 



(4.32) 



D, 



a n-l a n 



4 



Hq — S n o 



+ n(n+ 



(4.33) 



So, we should build IS for D n and _E ra and from them for H . Following the algebra, we 
find, for example for E n 

d 



E n = D n+ i + V n , 



V x = r — 1/\ 
n dr"' 



where 



e -r 2 r 4n+5 



(4.34) 



(4.35) 



n [A-2/ r e- 2 x 4 (™ +1 )dx]' 

Hence, if A < 0, then V*, v x and the respective eigenfunctions corresponding to the 
eigenvalue e n = (n + l)(n + 3/2), which become 



r d 
2d7 



(4.36) 



have not singularities and go to zero at infinity. Moreover, IS for H may be found from 
QOg ) using (|gD to give 



H-e n+h0 +e n + V n x , 



and writting the last equation as 



H n+lX 



n+l 



H + jV^ 



Xn+l — £n+l,0 Xn+l) 



(4.37) 



(4.38) 



where Xn+i = ipn, we note that H% is a family of semi-isospectral hamiltonians to H . 
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5. Concluding Remarks 

As we have shown, the existence of isospectral families of potentials is associated with 
the general factorization of the Hamiltonian H, even in the case when the operators 
which factorize to H are not ladder operators, and regardless H may has a continuos 
spectra {A n }. Also, there are IS when exist ladder operators, even though H is not 
explicitly factorized by them, but in this case, as we have shown above, they are not 
associated with H in a natural way. 

In the most general case, the eigenvectors {0 n } of the new hamiltonians 
do not inherit the property of satisfying the boundary conditions that the original 
eigenfunctions {</?„} satisfy, just as happened in the case of the free particle in one 
dimension, and when a^~^ = a^l\, {<f) n } will have the integrability properties of {<p n }- 
Moreover, the way to determinate the IS, is not unique, as in the case of the isotropic 
harmonic oscillator. In general, it is possible that exist more than one kind of IS 
associated to one operator, depending on how many independent ways to factored it 
can be found, and with that, the structures like the coherent states for the harmonic 
oscillator's isospectral family, related to the IS should be studied very carefully. 

On other hand, the algorithm discussed here has an strictly mathematical sense. 
Sometimes, as it happen in the case of the free particle, the new eigenvectors could 
have singularities and have not physical interest, but still this kind of functions could 
be important if we have bounded potentials, therefore they should not be rejected. 
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